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Abstract 

Recent theoretical and experimental investigations of coherent feedback control, 
the feedback control of a quantum system with another quantum system, has raised 
the important problem of how to synthesize a class of quantum systems, called the 
class of linear quantum stochastic systems, from basic quantum optical components 
and devices in a systematic way. The synthesis theory sought in this case can be 
naturally viewed as a quantum analogue of linear electrical network synthesis theory 
and as such has potential for applications beyond the realization of coherent feed- 
back controllers. In earlier work, Nurdin, James and Doherty have established that 
an arbitrary linear quantum stochastic system can be realized as a cascade connec- 
tion of simpler one degree of freedom quantum harmonic oscillators, together with 
a direct interaction Hamiltonian which is bilinear in the canonical operators of the 
oscillators. However, from an experimental perspective and based on current meth- 
ods and technologies, direct interaction Hamiltonians are challenging to implement 
for systems with more than just a few degrees of freedom. In order to facilitate more 
tractable physical realizations of these systems, this paper develops a new synthesis 
algorithm for linear quantum stochastic systems that relies solely on field-mediated 
interactions, including in implementation of the direct interaction Hamiltonian. Ex- 
plicit synthesis examples are provided to illustrate the realization of two degrees of 
freedom linear quantum stochastic systems using the new algorithm. 



1 Introduction 

Linear quantum stochastic systems (see, e.g., [H 121 El H]) arise as an important class of 
models in quantum optics [3] and in phenomenological models of quantum RLC circuits 
[B]. They are used, for instance, to model optical cavities driven by coherent laser sources 
and are of interest for applications in quantum information science. In particular, they 
have the potential as a platform for realization of entanglement networks [3 El E] and to 
function as sub-sytems of a continuous variable quantum information system (e.g., the 
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scheme of [10J ) . Recently, there has also been interest in the possibilities of control of 
linear quantum stochastic systems with a controller which is a quantum system of the 
same type [TTI |2"] |3| [T2] (thus involving no measurements of quantum signals), often referred 
to as "coherent-feedback control", and an experimental realization of a coherent-feedback 
control system for broadband disturbance attenuation has been successfully demonstrated 
by Mabuchi [12J. Linear quantum stochastic systems are a particularly attractive class 
of quantum systems to study for coherent control because of their simple structure and 
complete parametrization by a number of matrix parameters. 

A natural and important question that arose out of the studies on coherent control 
is how an arbitrary linear quantum stochastic system can be built or synthesized in a 
systematic way, in the quantum optical domain, from a bin of quantum optical components 
like beam splitters, phase shifters, optical cavities, squeezers, etc; this can be viewed as 
being analogous to the question in electrical network synthesis theory of how to synthesize 
linear analog circuits from basic electrical components like resistors, capacitors, inductors, 
op-amps, etc. The synthesis problem is not only of interest for quantum control, but is 
a timely subject in its own right given the current intense research efforts in quantum 
information science (see, e.g., [13]). These developments present significant opportunities 
for investigations of a network synthesis theory (e.g., [H] for linear electrical networks) 
in the quantum domain as a significant direction for future development of circuit and 
systems theory. In particular, such a quantum synthesis theory may be especially relevant 
for the theoretical foundations, development and design of future linear photonic integrated 
circuits. 

As a first step in addressing the quantum synthesis question, Nurdin, James and Do- 
herty [I] have shown that any linear quantum stochastic system can, in principle, be 
synthesized by a cascade of simpler one degree of freedom quantum harmonic oscillators 
together with a direct interaction Hamiltonian between the canonical operators of these 
oscillators. Then they also showed how these one degree of freedom harmonic oscillators 
and direct bilinear interaction Hamiltonian can be synthesized from various quantum op- 
tical components. However, from an experimental point of view, direct bilinear interaction 
Hamiltonians between independent harmonic oscillators are challenging to implement ex- 
perimentally with current technology for systems that have more than just a few degrees 
of freedom (possibly requiring some complex spatial arrangement and orientation of the 
oscillators) and therefore it becomes important to investigate approximate methods for 
implementing this kind of interaction. Here we propose such a method by exploiting the 
recent theory of quantum feedback networks [15]. In our scheme, the direct interaction 
Hamiltonians are approximately realized by appropriate field interconnections among os- 
cillators. The approximation is based on the assumption that the time delays required 
in establishing field interconnections are vanishingly small, which is typically the case in 
quantum optical systems where quantum fields propagate at the speed of light. 

The paper is organized as follows. Section [2] provides a brief overview of linear quantum 
stochastic systems and the concatenation and series product of such sytems. Section [3] 
recalls the notion of a model matrix and the concatenation of model matrices, while Section 
[4] recalls the notion of edges, their elimination, and reduced Markov models. Section [5] 
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reviews a prior synthesis result from |4J. Section [6] presents the main results of this paper 
and an example to illustrate their application to the realization of a two degrees of freedom 
linear quantum stochastic system. In Section [7] the special class of passive linear quantum 
stochastic systems is introduced and it is shown that such systems can be synthesized by 
using only passive sub-systems and components. Another synthesis example for a passive 
system is also presented. Finally, Section [8] offers the conclusions of this paper. In order 
to focus on the results, all proofs are collected together in the Appendix. 



2 Linear quantum stochastic systems 

This section serves to recall some notions and results on linear quantum stochastic systems 
that are pertinent for the present paper. A relatively detailed overview of linear quantum 
stochastic systems can be found in [I] and further discussions in [TJ, [2j [161 E] , thus they will 
not be repeated here. 

Throughout the paper we shall use the following notations: i = \J — 1, * will denote the 
adjoint of a linear operator as well as the conjugate of a complex number, if A = [a^] is 
a matrix of linear operators or complex numbers then A# = [a* k ], and is defined as 
A* = (A#) T , where T denotes matrix transposition. We also define, 3?{A} = (A + A&)/2 
and = (A — A#)/2i and denote the identity matrix by I whenever its size can be 

inferred from context and use I n to denote an n x n identity matrix. Similarly, denotes 
a matrix with zero entries whose dimensions can be determined from context, while Qmxn 
denotes a matrix with specified dimension m x n with zero entries. Other useful notations 
that we shall employ is diag(M 1; M 2 , . . . , M n ) (with M 1 , M 2 , . . . , M n square matrices) to 
denote a block diagonal matrix with Mi, M 2 , . . . , M n on its diagonal block, and diag n (M) 
(M a square matrix) denotes a block diagonal matrix with the matrix M appearing on its 
diagonal blocks n times. 

Let qi, pi, q 2 ,P2, Q n ,p n be the canonical position and momentum operators of a 
many degrees of freedom quantum harmonic oscillator satisfying the canonical commutation 
relations (CCR) [qj,Pk\ = 2ib^. The integer n will be referred to as the degrees of freedom 
of the oscillator. Letting x = (qi,Pi,q 2 ,P2, ■ ■ ■ ,Qn,Pn) T then these commutation relations 
can be written compactly as: 

xx T — (xx T ) T = 2iQ, 



with = diag n / 2 (J) and J 



1 
-1 



Here a linear quantum stochastic system G is 



a quantum system defined by three parameters: (i) A quadratic Hamiltonian H = \x T Rx 
with R = R T e M nxn , (ii) a coupling operator L = Kx where K is an m x n complex matrix 
and (iii) a unitary m x m scattering matrix S. We also assume that the system oscillator 
is in an initial state with density operator p. For shorthand, we write G = (S, L, H) or 
G = (S,Kx, \x T Rx). The time evolution, in the interaction picture, X(t) of x (t > 0) is 
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given by the quantum stochastic differential equation (QSDE): 



dX(t) = AX(t)dt + B 



dA(t) 
dA(t)# 

dY(t) = CX(t)dt + DdA(t), 



X(0) = x. 



(1) 



with 



A = 2Q{R + Q{K^K})- B = 2iQ[ -K^S K T S# ]; 
C = K; D = S, 

where A(t) = (Ai(t), . . . , A n (t)) T is a vector of input vacuum bosonic noise fields and 
Y(t) = (Yi(t), . . . ,Y n (t)) T is a vector of output fields that results from the interaction of 
A(t) with the harmonic oscillator. Note that the dynamics of X(t) and Y(t) are linear. We 
refer to (A, B, C, D) as the system matrices of G. For the case when n — 1, we shall often 
refer to the linear quantum stochastic system as a one degree of freedom (open quantum 
harmonic) oscillator. 

Elements of A(t) and Y(t) may be partitioned into blocks. For example, A(t) may be 
partitioned as A(t) = (A ri (t) T , . . . , A Tn . m (t) T ) and Y(t) as Y(t) = (Y S1 (t) T , . . . , F s „ out (t) T ) T , 
where A r .(t) and Y Sk (t) are vectors of bosonic input and output field operators of length 

m. We refer to n r . and n Sk 



n r . and n Sk , respectively, such that YlT=i n 



ETLout 
k=l 



J = L ■ -I j /L^lk=L ' ~ s k 

as the multiplicity of A r ,(t) and Y Sk (t), respectively. It is important to keep in mind that 
the sum of the multiplicities of all input and output partitions sum up to m, the total 
number of all input and output fields. With this partitioning, a linear quantum stochastic 
system G may be viewed as a quantum device having n- m input ports and n out output ports 
as illustrated in Figure [TJ The multiplicity of a port is then defined as the multiplicity of 
the input or output fields coming into or going out of that port. 




Figure 1: A linear quantum stochastic system viewed as a device with n- m input ports and 
n out output ports. 

Let us also recall the notion of the concatenation (EB) and series products (<) between 
open quantum systems developed in [TF]. The concatenation product G\ EB G2 between 
G\ = (Si, Li, Hi) and G2 = (S2, L 2 , H 2 ) defines another quantum linear system given by: 

G 1 EB G 2 = (diag(Si, S 2 ), [ L\ L\ ] T , H x + H 2 ), 

while if Gi and G 2 have the same number of input (and output) fields the series product 
G 2 < G± defines another linear quantum system given by: 

G 2 < G x = (S 2 Si, L 2 + S 2 Li, H 1 + H 2 + ^{LlS^}). 
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Note that in the definition of the concatenation and series product, it is not required that 
the elements of Li and L 2 and of Hi and H 2 are commuting. That is, G\ and G 2 may 
possibly be sub- components of the same system. Also both operations are associative, 
but neither operations are commutative. By the associativity, the operations EH" =1 Gj = 
G\ EH G 2 EH . . . EH G n and G n < . . . < G 2 < G\ are unambiguously defined. 

The concatenation product corresponds to collecting together the parameters of G\ 
and G 2 to form one larger concatenated system, and the series product is a mathematical 
abstraction of the physical operation of cascading G\ onto G 2 , that is, passing the output 
fields of Gi as the input fields to G 2 . The cascaded system is then another linear quantum 
stochastic system with parameters given by the series product formula. 



The model matrix and concatenation of model ma- 
trices 



The system G = (S, L, H) can also be represented using a so-called model matrix 
This representation will be particularly useful for the goals of the present paper. For G 
the model matrix representation M(G) is given by (the partitioned matrix): 



M(G) 



-iH - \lSL -tfS 



S 



(2) 



with the understanding that if L is partitioned as L = (Lj, L\, • • • , L£ ) t with Lj G rrij x 2 
and J2]j=i m j = m an d S is partitioned accordingly as S = [5 , jfc]j=i,..., nout ,A;=i,...,n iI1 with 
S jk e C m ' xm * and J2k=i 

m' k = m, then the model matrix above can be expressed with 
respect to this partition as: 



M(G) 



2 £uj=l "j^j £uj=l A 1 

L\ on 



E 



"out T t C 

j=l ^j D jn m 
Sln in 



J n out 



Sn 



(3) 



Also, with respect to a particular partitioning of M(G), it is convenient to attach a 
unique label to each row and column of the partition. For example, for the partitioning ^ 
we may give the labels sq,si,..., s„ out for the first, second, n out + 1-th row of M(G), 
respectively, and r ,ri, . . . ,r n . n for the first, second, n in + 1-th column of M(G), respec- 
tively. Moreover, with respect to this labelling (and analogously for any other labelling 
scheme chosen), elements of the blocks are denoted as: 



-j ^out '"OUt 

M soro (G) = -iH - £ 5>tL i5 M s r k {G) = -J^LtSjH, k > 0; 

j=i j=i 
M s . ro (G) = Lj, j > 0; M SjTk = S jk , j, k>0. 



"out 
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Since M(G) is another representation of a physical system described by G, and G can be 
identified directly from the entries of M(G), we will often omit the G and for brevity write 
a model matrix simply as M and denote its entries by M a p, with a ranging over row labels 
and (3 ranging over the column labels. Thus, we also refer to the triple S, L, H in Q as 
parameters of the model matrix M. 

Several model matrices can be concatenated to form a larger model matrix. Such a con- 
catenation corresponds to collecting together the model parameters of the individual ma- 
trices in a larger model matrix and is again denoted by the symbol EEL If G\ = (Si, L±, Hi) 
and G2 = (S2, L 2 , H 2 ) then the concatenation M(G\) EB M(G 2 ) is defined as: 



M(Gi) EB M(G 2 



-iHi - iH 2 



M(Gi ES G 2 ). 



1 ft T 
Li 
L 2 



l rt 7- 



-L\Si 
Si 




—L2S2 


s 2 



4 Edges, Elimination of Edges, and Reduced Markov 
Models 

Following [IS], a particular row partition labelled Sk with k > in a model matrix can 
be associated with an output port Sk (having multiplicity n Sk ) while a particular column 
partition rj with j > can be associated with an input port rj (having multiplicity n r ). 
In a system which is the concatenation of several sub-systems, it is possible to connect an 
output port Sk from one sub-system to an input port r 3 - of another sub-system (possibly 
the same sub-system to which the output port belongs) to form what is called an internal 
edge denoted by (sk,rj). For this connection to be possible, the ports Sk and rj must 
have the same multiplicity. Such an edge then represents a channel from port Sk to port 
rj. All ports which are connected to other ports to form an internal edge or channel are 
referred to as internal ports and fields coming into or leaving such ports are called internal 
fields. All other input and output ports that are not connected in this way are viewed as 
having semi-infinite edges (since they do not terminate at some input or output port, as 
appropriate) and are referred to as external ports and the associated semi-infinite edges 
are referred to as external edges. Fields coming into or leaving external ports are called 
external fields. From a point of view in line with circuit theory, one may think of a linear 
quantum stochastic system as being a "node" on a network and quantum fields as quantum 
"wires" that can connect different nodes. 

In any internal edge (sk,rj), there is a finite delay present due to the time which is 
required for the signal from port Sk to travel to port rj. As a consequence of these finite 
time delays, concatenated systems with internal edges cannot be represented by a Markov 
model such as presented in Section [2] (see [I] and the references therein for an overview of 
Markov models). However, as shown in [TSJ, the non-Markov model converges to a reduced 
Markov model in the limit that the time delay on all internal edges go to zero. That is, for 
negligibly small time delays, the reduced Markov model acts as an approximation of the 
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non-Markov model. In particular, such a reduced model serves as a powerful approximation 
of quantum optical networks in which signals travel at the speed of light and the time 
delay can be considered to be practically zero if the internal input and output ports are 
not extremely far apart. We recall the following results: 

Theorem 1 | I15L Theorem 12 and Lemma 16] Let T( Sk ,rj), j,k > 0, be the time de- 
lay for an internal edge (sk,rj) and assume that I — Skj is invertible. Then in the limit 
that Tr Sk ,r ) I 0, M(G) with the edge (sk,rj) connected reduces to a simplified model ma- 
trix M re( j with input ports labelled r Q ,ri, . . . ,rj_i,rj + i, . . . ,r„ out and output ports labelled 
s , si, . . . , Sk-i, Sfc+i, • • • , s„ in (i.e., the connected ports rj and Sk are removed from the la- 
belling and the associated row and column removed from M(G) ). The block entries of M re( j 
are given by: 

(M red ) Q/3 = M af3 + M arj (l - S^M^, 

with a e {s , Si, ... , s nout }\{sk} and (3 G {r ,r%, . . . , r nin }\{rj}. M Ted is the model matrix 
of a linear quantum stochastic system G re d with parameters: 

(Siedjpq = S p g + S p j(I Skj) Skq 

(Lredjp = L p + S p j(I — Skj) 1 Lk 
"out 

# red = H + J2 %{LtS pj (I - S k] )~ l Lp), 
P =i 

for all p e {1,2, ... , n ont }\{k} and q e {1,2,..., n in }\{j}. 

Several internal edges may be eliminated one at a time in a sequence leading to a 
corresponding sequence of reduced model matrices. The following result shows that such 
a procedure is unambiguous: 

Theorem 2 | I15L Lemma 17] The reduced model matrix obtained after eliminating all 
the internal edges in a set of internal edges one at a time is independent of the sequence 
in which these edges are eliminated. 

Suppose that M can be partitioned as: 



iH- l/2LfLi - l/2LtL e -L\S i{ - L\S ei -LtS^-L^ 



i 



L[ Sa S{ 



1C 



T q q 

-^e '-'ei ^ee 



(4) 



where the subscript i refers to "internal" and e to "external". That is, parameters with 
subscript i or ii pertain to internal ports, those with subscript e or ee pertain to external 
ports, while S[ e and pertain to scattering of internal fields to external fields and vice- 
versa, respectively. Interconnection among internal input and output ports can then be 
conveniently encoded by a so-called adjacency matrix. Let n\ denote the total multiplicity 
of internal input and output ports and let us view a port with multiplicity k as k distinct 
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ports of multiplicity 1. Suppose that these multiplicity 1 ports are numbered consecutively 
starting from 1, then an adjacency matrix 77 is an n\ x n\ square matrix whose entries are 
either 1 or with rj(j, k) = 1 (J, k e {1, 2, . . . , n;}) only if the j-th output port and the k-th 
input port form a channel or internal edge. Note that at most only a single element in any 
row or column of rj can take the value 1. Internal edges can be simultaneously eliminated 
as follows: 

Theorem 3 [15, Section 5] Suppose that M has a partitioning based on internal and 
external components as in and that connections between internal ports have been encoded 
in an adjacency matrix 77. If (77 — S 1 ;;)^ exists, then the reduced model matrix M re( j after 
simultaneous elimination of all internal edges has the parameters: 

5rcd = S ce + S e i(f] So) S{ e 

L TC d — L c + S e i(r) — Sii) 1 L[ 

H rcd = H+Y, ~ S^Li}. 

i=i>e 



5 Prior work 



Let G = (I m , Kx, \x T Rx) be a linear quantum stochastic system with K = [ K\ ... K n ] G 

C mx2n and R = jk , £ M 2nx2n ^ = ^T) Let Xj = (q^p-) T SO that X = 

(x± , . . . , x^) T and x satisfies the commutation relations of Section [2J Then the following 
result holds: 

Theorem 4 j^j Let Gk = (I,K k Xk, ^xlRkk^k)- Then 

G = (0, 0, H d ) ffl (G n < . . . < G 2 < d), 
where H d = ^ Y2 =j+ i K ™d E% = x T k (Rj k - 



Aft) 



G = (I,L,H) 



Hf 2 



Hi 



£ 



*-o G, o 



o G 7 o 



G„ o 



Hi 



Y(t) 



Figure 2: Realization of G = (J, L, H) via a cascade connection of one degree of freedom 
quantum harmonic oscillators Gj and bilinear direct interaction Hamiltonians Hjk between 
every pair of oscillators 

The theorem says that G can in principle be constructed as the cascade connection G n < 
. . . < G2 < G\ of the one degree of freedom oscillators Gi, G2, ■ ■ ■ , G n , together with a direct 
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G 

~G~ 

A r Jt) 
A r,ft> r 



-o r u G„ s " o- 

-O'n G 12 s '-'0 



Ys.fi) 



G, 

■A r Jt) 

A rJt) 



~O r n G 2I S2 '°- 

-0*22 G„'«o- 



Y s Jt) I 



G„ 

Aft) r 



-0>;„ G„, s -0- 



Y s ft) I 



Y s ft>: 



• • • 



4-. (0 



~On„ G,„ J '»0 



ft) 



_i 



A ft; 



v^ 1 



^ ft). 



»o- 



Y s (t) 



Figure 3: The concatenation decomposition of G as G 



ffl* 



En 
, fc= 



Go 



ifc- 



bilinear interaction Hamiltonian H d which is the sum of bilinear interaction Hamiltonians 
Hj k between each pair of oscillator Gj and Gt- This construction is depicted in Figure [2 
It was then shown how each Gk can be constructed from certain basic quantum optica 
components and how H^ k can be implemented between any pair of oscillators. However, 
a drawback of this approach, based on what is feasible with current technology, is the 
challenging nature of implementing H d . This may possibly involve complex positioning 
and orientation of the oscillators and thus practically challenging for systems with more 
than just a few degrees of freedom. Although advances in experimental methods and 
emergence of new technologies may eventually alleviate this difficulty, it is naturally of 
immediate interest to also explore alternative methods of implementing this interaction 
Hamiltonian, at least approximately. In the next section, we propose such an alternative 
synthesis by exploiting the theory of quantum feedback networks that has been elaborated 
upon in preceding sections of the paper. 



6 Main synthesis results 

For j = 1, . . . , n, let G jk = (S jk , L jk , 0) for k = 1, . . . , n, k ^ j, and Gjj = (%, Lj, Hj), 
with Sjk G C C: > kXCjk , Cjk = Ckj and Cjj = m, Lj k = Kj k Xj, and Hj = ^xjRjXj with 
Rj = Rj G M 2x2 . Here Xj = (qj,Pj) T is as defined in the previous section. Let Gj = 
S k=1 Gjk for j = 1, . . . , n, and note that Gj = (Sj, Lj, Hj) with Sj = diag(5^i, Sj2, ■ ■ ■ , Sj n ), 
Lj = (Lj 1 ,Lj 2 , • • • , Lj n ) T , and Hj as already defined. 

Consider now the model matrix M for the concatenated system EB" =1 Gj, see Figure [3j 
With respect to the natural partitioning of M induced by the Gj k s (via their concatena- 
tion), we label the first n + 1 rows of M as s o, s n> • • • > s in, the nex t n rows as s 2 i, • • • , s 2n , 
the n rows after that s$i, S32, S33, . . . , S3„ and so on until the last n rows are labelled 
s n i, s n 2, . . . , s nn . Similarly, we label the first n + 1 columns of M as r Q0 , ru, . . . , ri n , the 
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next n columns as T2i, ■ ■ ■ , ?"2 n , the n columns after that r$i, . . . , r% n and so on until the 
last n columns r nl ,r n2 • • • ,r nn . On occasions, we will need to write a bracket around one 
of both of the subscripts of r or s (e.g., as in S( n _i)fc or r^+i)(k+i))- 



Theorem 5 Let the output port Sjk be connected to the input port r^ to form an internal 

-Sjk I 

is invertible Vj, k — 1, . . . , n, j 7^ k, then the reduced model matrix M re d obtained by allowing 
the delays in all internal edges {e^, j 7^ k} go to zero has parameters given by: 



edge/channel = (sjk, r k j) for all j, k = 1, . . . , n, j 7^ k. Assuming that 



5red 
-^red 



diag(S , n, S22, • • • , S nn ) 

n— 1 n 



Hi 



red 



E^+E E 

j=i k=j+i 



k=l 





-1 


L jk 




—Skj 1 




. Lkj 


} 



The proof of the theorem is given in Appendix |A} We will also exploit the following 
lemma: 



Lemma 6 For any real 2x2 matrix R and unitary complex numbers S12 and S21 satisfying 
S12S21 7^ 1, there exist 1x2 complex matrices K\ and K 2 such that R — 1 _g 12 g K1K2 + 
^f" KfK*} = 0. In fact, a pair K\ and K 2 satisfying this is given by K\ — [ k in ] 
with k an arbitrary non-zero real number and K2 = 2i[ 1 ][ —KlA* K^ A ]~ l R, where 
A = 2 | 1 _g i2 g^| 2 . Or, alternatively, K 2 = [ k in ] and 
K x = 2i[ 1 ][ KlA -KjA* Y l R T . 



See Appendix [B] for a proof of the lemma. As a consequence of the above theorem and 
lemma, we have the following result: 

Corollary 7 Let Cjk = Ckj = 1 whenever j 7^ k, and Cjj = m for all j, k = 1, . . . , n. Also, 



let Si 



jj 



Im, S k j = e l&k > and S jk 



e Jfe with 9kj,0jk G [0, 2n) satisfying 6kj + Ojk 7^ 7 
Kjj = Kj andRj = Rjj - 2sym fefc=i,fc^ ^{ is] k s kj K ]k K jh}^ (Kj and R 3J given), where 
sym(A) = |(v4 + A T ), and the pair {Kjk, Kkj) (j 7^ k) be given by: 

Kjk 
K k j 



[ Kjk IKjk \ 

2i[l 0)[-K} k A* k Kj k A jk )-\R jk -${KjK*}), 



(5) 



or 



K 



kj 



Kjk IKjk 



K 



where Rjk 



jk 



R kj G 



2i[l O^K^Ajk -K^A* k ]- 1 (R 1 k-^{KjK*}) T , 



p2x2 



'kj 

Skj 



(6) 



jk 



\l—SkjSjk\^ 



and Kj k is an arbitrary non-zero real constant 



for all j, k. Then the reduced Markov model G re d = (S ie d, L ve d, H ie d) has the decomposition 
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Cred = fflfc =0 G r ed,k with G rc d,o = (0, 0, # rcd ) and G rc d,k = (S kk ,L kk ,0) for k = l,...,n. 
Moreover, the network G net = (S net , L net , H net ) formed by forming the series product of 
Gmd,n < • • • < C ro d,2 < G'red.i within the concatenated system G re d and defined by G net = 
G re d,o ffl {G re d,n < • • • < G>ed,2 < C re d.i) a linear quantum stochastic system with parameters 
given by: 



Snet 
L ve d 



Kx, K =[ K x K 2 



K n 



red 



~-x T Rx, R = [R 



■jk\j,k=l, 



In other words, G nct realizes a linear quantum stochastic system with the above parameters. 



Remark 8 Note that the series connection G re d,n < • • • < G rc d i2 < G Te d 7 i can be viewed as 
forming and futher eliminating the internal edges {(sj-k, r( k +i)(k+i))) k = 1, . . . , n — 1} in 
G Te d or G and, hence, is in essence a special case of a reduced Markov model [Th^ . 

The proof of the corollary is given in Appendix [Cj The corollary shows that an arbitrary 
linear quantum stochastic system (J, Kx, ^x T Rx) can be realized by a quantum network 
G net constructed according to Theorem [5] and the corollary, with an appropriate choice of 
the parameters Rj, Kj k and S k j (j ^ k). From here, any system (S, Kx, \x T Rx) can then 
be easily obtained as (S, Kx, \x T Rx) = (J, Kx, t,x t Rx)<(S, 0, 0). That is, as a cascade of a 
static network that realizes the unitary scattering matrix S and the system (J, Kx, \x J Rx) 

(see (mm). 

Example 9 Consider the realization problem of a two degrees of freedom linear quantum 
stochastic system G sys = (I,Kx, l/2x T Rx) (x = (qi,pi, q2,P2) T ) with 



K=[K X K 2 

Rn 



[ 3/2 i/2 1 % ) (Kj G C lx2 , j = 1,2 



R 



R 12 



R\2 
R22 



2 

0.5 
1 
1 



0.5 
3 

-1 
-1 



1 -1 

1 



1 





(Rjk £ 



D 2x2 



, j,Ar = l,2) 



Let xi = {qi,pi) T and x 2 = (q2,P2) T - Define G\ = ( diag(J m , Si 2 ), 



K 12 



1 5 9*^1 



and G 2 = diag(S , 2 i, I mJ . 



mined. Choose 612 = 



K21 
K22 
and 821 
Then set S 12 = e lf)l2 = I, S 21 
[ 1 i ]. Compute Ai 2 = 

we set K\2 = [ K42 ^12 
1 0.5 
0.5 2 



X2, \x 2 R 2 x 2 y with certain parameters still to be deter- 
/2, so that 9\ 2 + 621 7^ as required in Corollary^ 



71 



2 S21—SI2 



Ri 



and R 2 = — 



1, K n = Kx = [ 3/2 1/2 ] and K 22 = K 2 = 
[_ S|lS ' 21 p = —1 + i and set Kx 2 = 1. Then by Corollary^ 
= [ 1 i] and compute K 21 = [ 1.25 - 0.25z 1.75 + 0.75z ], 

1 Thus, we have determined all the param- 



2 2.625 

eters of G± and G 2 . Labelling the ports of G\ and G2 according to the convention adopted 
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in this section, G sys can be implemented by concatenating G\ and G 2 and eliminating the 
internal edges (si 2 ,r 21 ) and (s 2 i,r 12 ) to form G re( j and then eliminating the edge (sn,r 22 ) 
(cf. Remark^ to obtain G Qet as an approximate realization of G sys . This realization is 
illustrated in Figure [^| G\ and G 2 can then be physical realized in the quantum optics 
domain following the constructions proposed in ^j. Using the schematic symbols of FJ]/ 7 a 
quantum optical circuit that is a physical realization of G sys is depicted in Figure [5| 




Figure 4: Realization of G sys via a quantum feedback network. In the top figure the internal 
edges (si 2 , r 2 i) and (s 2 i, r i2 ) are formed and eliminated to obtain a reduced Markov model 
Gred- Then in the bottom figure a series (cascade) connection is formed by the eliminating 
the internal edge (%, r 22 ) to realize G sys 



7 Synthesis of passive systems 



Let us now consider a special class of linear quantum stochastic systems that we shall refer 
to as passive linear quantum stochastic systems, for a reason that is explained below, and 
show that any such system can be built from passive components. This class of systems 
has also been considered in, e.g., 

For k = 
a = (ai, . . . 



. . . , n let a*; = ik+^Pk k e annihilation operators for mode k and define 
) T . Then a satisfies the CCR 



Moreover, note that: 



where 



[ 



a' a 



T 



] 



diag(7„, -I n ). 



a 
a 



# 



a 






_ a* _ 






i\ 





( 





1 
2 


i\ ( 



,r; 



[ 












12 



Am 



Auxiliary X 



Figure 5: A quantum optical circuit that realizes G sys according to the quantum feedback 
network in Figure 4 the schematic symbols are taken from jl]. Here aj = gj is cavity 
mode of the optical cavity around which the physical realization of the oscillator Gj is 
1,2. The dashed box labelled Gj is the part of the circuit realizing Gj, and 



based, 



J 



the input and output ports of Gj are indicated by their respective labels. The value of the 
parameters of devices in this circuit can be determined according to jl]. 



We also make note that 



E 

E# 



-i -l 



x 



E 

E# 



2 [ £t E T ] and therefore we have 
= 2 [ £t s T ] 



An n degrees of freedom system G = (S, Kx, \x Rx) is said to be passive if we may 
\a) Ra + c and L = i^a; = Ka for some real constant c (arising due 



\x T Rx 



write 

the commutation relations among the canonical operators), some complex n x n Hermitian 
matrix R and some complex mxn (m here again denotes the number of input and output 
fields in and out of G) matrix K. Here the term passive for such systems is physically 
motivated. By the given definition, the Hamiltonian H contains no terms of the form c\Oj, 
c 2 a* k 2 , c 3 aj(ik and c^a*a\ and the coupling operator L contains no terms of the form c 5 a* k 
(here c\, . . . , C5 denote arbitrary complex constants), with the indexes j and k running over 
1, . . . ,n. These terms are precisely the terms that require an external source of quanta 
(i.e., an external pump beam) to realize (see, e.g., Jl]) and cannot be implemented only 
by passive components like mirrors, beamsplitters and phase shifters. The absence of such 
"active" terms is the reason we refer to this class of systems as passive. 
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Let us write \a) Ra and Ka in the following way: 



a) Ra 



at a T ] 







2-R nX n 
nxn 2 



a 



1 n 



1 1 



if a = K [ I ] 



a 



Therefore, H = \x T Rx with R = ^{Ti^RT,}, and L = Kx with = KT,. It is easy to see 
from this that for any passive system the block diagonal elements Rjj must be diagonal 
matrices of the form Xjl 2 with Xj G K, the off diagonal block elements Rjk are real 2x2 

C^j'fc Pjk 



matrices of the form R^ 



for some real numbers and (3^ and the 



coupling matrix Kj to Xj is of the form Kj = [ "fj £7^ ] for some complex number jj. One 
would hope that a passive system can be synthesized using purely passive components and 
the next theorem states that this is indeed the case. 

Theorem 10 Let G sys = (S,L,H) be passive. Then the systems {Gj;j = l,...,n} con- 
structed according to Corollary [7| are all also passive. 



The proof of the theorem can be found in Appendix [Dj We now conclude this section 
with an example of a passive system synthesis. 



Example 11 Consider the realization problem of a passive two degrees of freedom linear 
quantum stochastic system G sys = {I,Kx, l/2x T Rx) (x = (qi,pi, a 2 - l V'i) T ) with 



'sys 

K —[ Kx K 2 ] = [ -3 + i 

~ 2 

-Rn R12 2 



1 — 3i 1 % ] {Kj e C lx2 , j = 1,2); 



R 



R l2 R22 



1 4 
-4 1 
1-410 
4 10 1 



{Rjk e 



)2x2 



j,k = 1,2). 



Setting k 12 = I, S\ 2 = 1, 5*21 = i and K 12 = [ ^12 ««i2 ] = [ 1 i }, we obtain from 
CorollaryQ K n = K x , K 22 = K 2 , K 21 = [ 0.5 - 0.5i 0.5 + 0.5z ], R x = 2x2 and R 2 = 
0.5/2- An entirely passive optical circuit that realizes G sys is shown in Figure^ 



8 Conclusions 

In this paper we have developed a new method or algorithm for systematically synthesizing 
an arbitrary linear quantum stochastic system via a suitable quantum feedback network. In 
the synthesis, all interactions between the oscillators constituting the network are facilitated 
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Figure 6: A quantum optical circuit that realizes the passive system G sys of this example. 
The circuit consists of only passive optical components: mirrors and phase shifters. Here 
#11 = 7r — arctan(l/3) and 6*21 = — f , and the values of all other parameters of devices in 
this circuit can be determined according to [I]. 



only by quantum fields. In particular, it offers an alternative to the challenging task of 
implementing direct bilinear interaction Hamiltonians that was required in the approach 
of [1]. Moreover, it is shown that if the system is passive then with the new algorithm it 
can be realized using only passive components. It is clearly also possible to combine the 
present algorithm with that of [I] to form a hybrid synthesis method. 

Current interest in quantum information science, and quantum control in particular, 
provides much impetus for extending network synthesis theory to the quantum domain 
as a significant direction in which to further develop circuit and systems theory. The 
quantum synthesis results presented here could be of particular relevance for the design 
and systematic physical realization of photonics based monolithic linear open quantum 
circuits which may form an important sub-system of future continuous variable quantum 
information processing systems. 



A Proof of Theorem 5 



For this proof, it will be convenient to interchange some rows and columns of the model 
matrix M to form another model matrix M to avoid complicated book keeping and thus re- 
duce unnecessary clutter. Here by "rows" and "columns" we mean, respectively, block rows 
and block columns of M formed with respect to its specifed partitioning. This interchange 
is as follows. 

First, we permute rows of M such that the first 2n — 1 rows from top to bottom are the 
rows labelled (while column labels are kept fixed as they are) soo, Si2j S2i> Si3> S3i> • • • , sin> s ni, 
the next 2(n — 2) rows respectively are the rows labelled S23, S32, S24, S42, . . . , S2 n , s n 2, the 
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2(n — 3) rows after that are respectively the rows labelled S34, S43, S35, S53, . . . , ss n , s„3 
and so on in the same pattern until we get to the last n rows that are respectively 
those rows labelled Sn, s 22 , • • • , s nn . Call the intermediate matrix resulting from this row 
permutation M. Then fixing the row labels of M, we permute its columns such that 
the first In — 1 columns from left to right are respectively the columns of M labelled 
r oo, r i2, ^21, ^13, ^31, ■ • ■ , T\ n , r ni, the next 2(n — 2) columns are respectively the columns 
labelled r 2 3, r 32 , r 24 , r 42 , . . . , r 2n , r n2 , the 2(n — 3) rows after are respective the columns la- 
belled r 34 ; r 43; r 35 ; r 53 ; • • • j r 3n ; r n3 ; and so on in the same pattern until the final n columns 
that are respectively the columns labelled r 11; r 22 , . . . , r nn . The resulting matrix after this 
permutation of columns is M. 

It is important to note here that since the same permutation is applied to the rows 
and columns, M and M are model matrix representations of the same physical system, 
save for a mere rearrangement of the ordering or indexing of the fields and ports. That 
is to say that if M is the model matrix of G = (S, L, H) then M is the model matrix 
of G — (PSP T , PL, H) for some suitable constant real permutation matrix P, while it is 
clear that G and G are representations of the same physical system. Thus with the same 
internal connections made, a reduced model matrix for M is also a reduced model matrix 
for M, up to a possible relabelling of uneliminated ports. 

Let L = PL and S = PSP T . Then L can be partitioned as L = (Lf, L 7 ) T , where L\ 
is the first n(n — 1) + 1 rows of L, while L e is the last n rows of L. They are of the form: 

T t t 



Li 



T jT jT jT 

12) -^215-^13' -^31' 



T T jT jT jT jT jT 

^Ini n nli -^23' -^32' -^24' -^42' 



(L 

T 7 I 1 I 1 I 1 I 

-^(n-2)(n-l)) - tv (n-l)(n-2)) ^(n-2)ni - tj n(n-2)) ^(n-l)n> 



T 1 T 1 



T T T jT \T 

n n(n-l)) 



(Lji,L 



T 
22' 



T \T 



Similarly, 5* can be partitioned as 5* 



q.. q. 

q ■ q 



with 5*ii and S ee being block diagonal: 



Sn 



diag(S"i 2 , S21, 5*i3, 531, . . . , Si n , S n i, S23, S32, S24, S42, ■ ■ ■ , SWi, 5^2, • • • , 

S(n-2)(n-l)j 5( n _i)( n _ 

diag(5n, S22, ■■■,S r) 



Q Q Q Q Q Q 

>->(n-2)(n-l)> >->(n-l)(n-2)j >->(n-2)rc> >->n(n-2)> >~>(n-l)n, ^n(n-l) 



and S e i and Si e both being zero matrices. Then M has the partitioning of the form (|4j) by 
identifying Sn, Si e , S e i, S ee , L\ and L e with Sn, Si e , S e {, S ee and L; and L e , respectively. 
The reduced model matrix resulting from the subsequent simultaneous elimination of all 
internal edges (sjk, r^j) j, k — 1, . . . ,n, j ^ k, can be conveniently determined by using the 
adjacency matrix n defined by: 



V 



diag 



( 


J C12 ~ 




I C13 ' 




_h 21 




_h 31 










J C23 " 




J C24 " 


Ic , 




. Ic 32 




. h 42 




1 













* * 1 


Ic n2 





, . . . , 


Jo 




n — l)n 
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(Recall that Cjj. = Ckj). Hence, according to Theorem [3J the reduced model matrix M reo \ 



obtained after elimination of the internal edges {(sjk,r k j); j,k = 1, . 
parameters given by (recalling that S ci and S ie are zero matrices): 



*Sred 



red 



See + Sei(V Sii) S{ c — S ec , 

L e + S e i(f] — Sin) Li = L e 

n 

J2 H ^ + Y1 ^Aiiv - Sur 1 Li} 

k=l i=i,e 
n 

k=i 

n 

n n—1 n 

E*+EE ! 

j=i k=j+i 



fc=i 



L 



jk kj 



I 

—Skj 



-S jk 

I 



n,j 7^ k}, has 



-1 








. Lkj 


} 



Since M and M are model matrix representations of the same physical system and the 
external fields have the same ordering and labelling in both representations, the reduced 
model matrix of M re d and M re d of M and M, respectively, after elimination of internal 
edges (sjk,r k j), coincide. Hence, also the linear quantum stochastic systems G re d and G Tec \ 
associated with M and M, respectively, coincide. This completes the proof. □ 



B Proof of Lemma 6 



We begin by noting that 
S12 



K\K 2 + 



12*J21 



521 KfK*} = ^[-KlA* KfA 



1 — S12S21 



K, 



# 



with A = jz§^ ~ i-sts* 12 = 2 |i-4iS| 2 ( ex P loitin § the fact that SuSu = 1 = ftiSaii)- 
Now set K\ = [ k in ] for an arbitrary non-zero real constant K, and note that S21S12 7^ 1 
implies that A^O and: 



-K\A* K{ A 



-kA* kA 
ikA* ikA 



-1 -1 



1 



2ik 2 \A\ 2 



inA —kA 
-ikA* -kA* 



and therefore for any real matrix V, 2i[ —K\A* K^A } l V 



Z 
Z* 



for some complex 



row vector Z. Therefore, given any R we see that we may solve the equation 



[ -k\a* KfA 



K. 



# 



2iR, 
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for K2 and this solution is as given in the statement of the corollary. 

Alternatively, we could also have started by setting K 2 = [ n in } and analogously 
solving for K\ for a given R. It is then an easy exercise that the solution for K\ in this 
case is as stated in the corollary. □ 



C Proof of Corollary [7] 

With Cjfc, Sjk, Rjk and Kjk, defined in the statement of the corollary, 

from Theorem |5] and Lemma [6] we have that S re d = I n m, A-cd = {LJijL^, ■ ■ ■ ,L^ n ) T with 
Ljj — — fcCjOCj^ cincl 



n— 1 n 



j=l k=j+l 



3=1 

Expanding, we have: 



1 


-S jk 


-1 


Ljk 




1 







red 



■ikSkj 



2 E4^-E E Hr^ 

j=i j=i k=j+i ^ 3 

(Ljk-Ljk + SjkL^ jk L k j + SkjL^Ljk + Lj^Lfcj) j 
;K+2sym( ^ 



n— 1 n 



k=\.k£j 



SjkSkj 



fc— ^jk^k] 1 ~ ' '/A-' '/.-/ 



n— 1 n 



jE^ Ti? ^ +E E *J(%-9{JEpf?}) 



jk^kj 
X k 



3=1 



j=l fc=j+l 



n— 1 n 



Vi?* -EE x Pi K I K t 



j=l k=j+l 



Rj k . From this it is 



where sym(A) = + v4 T ), and R = [Rjk]j,k=i,...,n and Rkj 

clear that using the concatenation product we can decompose G re d as G ro d = (0, 0, i? re d) ffl 
ffl" =1 (7 m , L#, 0). Let G redi0 = (0, 0, H vcd ) and G rcdii = (J m , 0), j = 1, . . . , n. Now, using 
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the series product rule (cf. Section [2]), we easily compute that 

G ne t = Gredfi EB {G Ie d,n < ■ ■ ■ < G re d,2 < G re d,l) 

n— 1 n 



II — _L II \ 

j=i k=j+i ' 

(n— 1 n \ 
I m ,[Ki K 2 ... K n }x,J2 Yl X J^{ K I K *M 
i=i k=j+i ' 



3=1 k=j- 

I m , [ K i K 2 ■■■ K n ]x, ^x T Rx ) . 

Therefore, G net realizes a linear quantum stochastic system with parameters S nct , L net 
and H net , as claimed. □ 



D Proof of Theorem 10 



for some a'- fc , /3'- fc G E, whenever j ^ k. Let Kjk, Sj k and S k j be chosen 



By the passivity of G sys , we immediately see that K 33 = Kj is already of the form required 
for passivity of Gj, and the matrix Rj k — Q{KjKf} is a 2 x 2 real matrix of the form 

—f3j k aj k 

according to Corollary [7j Set Kj k according to the first equality of Then some 

straightforward algebra shows that K k j given by the second equality of is of the form 
K k j = [ 7^- ij' k j ] for some jL G C, just like K jk . The same holds true if one chooses 
the alternative of setting K k j according to the first equality of ^ and computing Kj k 
according to the second equality of ([6]). Finally, given this special form of Kj k , it is 
easily inspected that sym(S ; { 1 _ s 1 ^ s ^ Kj k Kj k }) is a diagonal matrix of the form Xjkh for 

some Xjk G K for all j, k,j ^ k (recall that sym(A) = |(^4 + A T )), and since Rjj is 
also diagonal of this form (again from the passivity of G sys ) we have that R 3 = Rjj — 
2 J2k=i k^j sym(^{ 1 _ s 1 fcSfc Kj k Kjk}) is again of the same form, for all j. Therefore, it now 
follows that each sub-system Gj of Section [6] with parameters determined according to 
Corollary [7] is passive. This completes the proof. 
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